arXiv:1502.07740vl [astro-ph.GA] 26 Feb 2015 


Mon. Not. R. Astron. Soc. 000, [I][8] (2002) Printed 27 February 2015 (MN IAT^X style file v2.2) 


The Contribution of Faint Galaxy Wings to 
Source-subtracted Near-infrared Background Fluctuations 

R.L. Donnerstein 1 * 

1 The University of Arizona, Tucson, AZ 85724, USA 


Accepted 2015 February 16. Received 2015 February 5; in original form 2014 July 11 


ABSTRACT 

The source-subtracted, 1.1 and 1.6 (im NICMOS images used in earlier analyses of the 
near-infrared Hubble Ultra Deep Field contained residual flux in extended wings of 
identified sources that contributed an unknown amount to fluctuation power. When 
compared to the original results, a reanalysis after subtracting this residual flux shows 
that mean-square and rms fluctuations decrease a maximum of 52 and 31 per cent at 
1.6 /im and 50 and 30 per cent at 1.1 /iin. However, total mean-square fluctuations 
above 0.5 arcsec only decrease 6.5 and 1.4 per cent at 1.6 and 1.1 /zir i, respectively. 
These changes would not affect any published conclusions based on the prior analyses. 
These results exclude previous suggestions that extended wings of detected galaxies 
may be a major contributor to the source-subtracted near-infrared background and 
confirm that most fluctuation power in these images must be explained by other means. 

Key words: infrared: diffuse background - methods: data analysis - cosmology: 
observations. 


1 INTRODUCTION 


Multiple theories have been offered to explain residual 
fluctuations in the near-infrared background (NI RB). Af- 
ter analysi ng images from Spit zer /IRAC ( Kashlinsk y et al.l 
120051 . 12 OI 2 I ') and AKARI /IRC dMatsumoto et al.ll201 il l . the 
authors concluded that NIRB fluctuations are unlikely to 
result from known galaxy populations and are consistent 
with the earliest star formation era. On the other hand, 
evaluation of NICMOS images from the Hubble Ultra Deep 
Field (NUDF) suggested that source-subtracted fluctua¬ 
tions result from normal galaxies at 2 < 8, with most 
power in the redshift range of 0.5-1.5 llThompson et all 
l200TaHbh . More recently, other explanations have been pro¬ 
posed to account for the majority of near-infrared fluctu¬ 
ations, including intrahalo light (IHL) from d ark-matter 
haloes dCoorav et al.ll2012bl : IZemcov et al.ll2014l 'l and early, 
di rect c oll apse black holes (DCBHs) as suggested by 
lYue et ah! (l2013blf . Conclusions are frequently based, in 
part, by the observation that NIRB fluctuations measured 
by different instruments over a broad range of wavelengths 
are consistently m uch high er t h an t hose p redicted for 
known populations dHelgason, Ricotti, fe Kashlinskvll2()l2l) 


or high-redshift galaxies dCoorav et al.l 12012a : Yue et all 


l2013al : iFernandez et al.l [2012). Fluctuation measurements 
considered in these studies include those derive d f rom 1.6 
(im NUDF observations as described by [Thompson et all 


J2007all . All of these interpretations depend upon the am¬ 
plitudes and shapes of observed fluctuations, which, in turn, 
are r elated to the completeness of the source- extraction pro¬ 
cess. iHelgason. Ricotti. fe Kashlinskvj (j2012i l appropriately 
noted that the source-sub tracted NUDF images analysed by 
iThompson et all d2007al l included residual flux in the ex¬ 
tended outer edges of galaxies (wings). This residual flux 
would impact fluctuation estimates to an unknown extent 
and could potentially affe ct conclusions der i ved fr om their 
spectra. As described by [Thompson et al] (|2006l ). source 
subtraction for these images was deliberatively conserva¬ 
tive to minimise the possibility of false sources appearing 
in their catalogue. This resulted in incomplete subtraction 
of detected objects with visible flux remaining in the wings 
of identified sources (Fig. [ljc). If these unsubtracted wings 
are a substantial component of measured fluctuations, con¬ 
clusions drawn from their result s migh t have to be recon- 
sidered. While iKashlinskv et all (l2005l l and lArendt et al.l 
d20ld) showed that residual wings contributed only small 
amounts to the power spectra of Spitzer/IRAC images, their 
methodology was limited by the resolution of this system. 
The higher resolution of i/u&We/NICMOS observations al¬ 
lows a more robust approach for assessing t hese effects. Thi s 
paper r eanalyses the NUDF images used bv lThompson et al.l 
(l2007arl to quantify the contributions of these unsubtracted, 
extended galaxy wings to measured fluctuations. 
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Figure 1. Image examples, (a): Original, full F160W image, 
(b): Source-subtracted, full F160W image. White box outlines 
enlarged region shown in the lower panels and in Fig. [3] (c): En¬ 
larged F160W image. Residual flux is evident in source wings. 
Black box outlines the source used for the detailed example shown 
in Fig. [2] (d): Enlarged F110W image. Source wings are still seen, 
but are not as prominent as in the F160W image. Panels (b)—(d) 
have the same linear stretch. Unless otherwise noted, masked pix¬ 
els in figures are set to black for illustration purposes. These pixels 
are assigned values of zero for fluctuation calculations. 


2 OBSERVATIONS AND DATA REDUCTION 

All images and tables used in this study are available onlineQ 
Detailed descri ptions of the proc e ssing of NUDF images are 
documented in iThompson et ah 12005 1. Only overviews of 
portions relevant to this study are presented. 


2.1 Image acquisition and processing 

NUDF observations cover a 144 x 144 arcsec area and are 
obtained with two filters centred on wavelengths of 1.1 
(F110W) and 1.6 (F160W) (jm. Field position on the camera 
is dithered by more than the average source size during sepa¬ 
rate acquisitions whic h results in a median ima ge dominated 
by zodiacal emission (IThompson et al1l2007al l. Subtracting 
this median image from the original images effectively re¬ 
moves smooth background components, including the zodi¬ 
acal light. Images are then combined in a drizzle procedure 
such that the original 0.2 arcsec pixels are converted to 0.09 
arcsec to match a 3 x 3 binning of the 0.03 arcsec Hubble 
Advanced Camera for Surveys (ACS) images. 


2.2 Source extraction 

The source e x tracti on procedure is described in 
IThompson et all (1200614 and is briefly summarised here. 


1 http: //archive.stsci.edu/pub/hlsp/udf/nicmos-treasury/ 


Using a process described by ISzalav. Connolly, fc Szokolvl 
l|l999l) . all four ACS and the two NICMOS bands are used 
to identify those pixels with a sufficient signal to noise ratio 
to be considered part of a real source. Source extraction 
is then done with SExtractor (SE) i n the t wo-image mode 
as described by iBertin fc Arnoutsl (Il996l 4. This process 
uniquely identifies individual sources and their locations 
on the image@ The ACS optical images have much longer 
integration times than NICMOS images and, therefore, 
identify sources not apparent in the NUDF images alone. 
This results in more complete source detection than would 
be possible using only NICMOS observations. A total of 
4702 objects are listed in the NUDF catalogue generated 
by SE0,' 4276 of which remain after trimming. Sources 
identified by this method are then subtracted (masked) 
from the image by setting their corresponding pixels to 
zero. As shown in Fig. [l] this process may fail to subtract 
outer portions of galaxies more prominent in the NICMOS 
images. Detector quantum efficiency is less at 1.1 /im than 
at 1.6 /im causing signals from the outer regions of galaxies 
to blend more quickly into background noise. This makes 
source subtraction for the F110W image in Fig. [ljl to 
falsely appear more complete than that for the F160W 
image (Fig. QJ). 


2.3 Fluctuation analysis 


Fluctuation an alyses are performed as d escribed in detail in 
Appendix A of IThompson et al. rt2007ai) and are similar to 


those used by others ( Kashlinskv et all 20051 : ICoorav et al.l 
2007). The original image in units of ADU s -1 is converted 
to Jy pixel -1 using standard procedures for NICMOS im¬ 
ages. Pixel values are then converted to W m -2 sr -1 by 
multiplying by the corresponding frequency for 1.1 or 1.6 
|im and dividing by the solid angle subtended by each pixel. 
Each source-subtracted image is then multiplied by a nor¬ 
malised weighting image to account for different pixel ob¬ 
serving times. Prior to Fourier transformation, the mean 
value of the unmasked pixels are subtracted from this re¬ 
sult to produce an image with a mean of zero. The two- 
dimensional Fourier transform of the resultant image, f(q), 
is computed using an FFT, where q is the wave vector. 
The value of the wavenumber, q, and the amplitude of the 
transform, |/(g)|, is calculated for each point in the re¬ 
turned array. A log-spaced vector of q values is generated 
and the average value of \f(q)\ 2 , P 2 , is computed for each 
of the bins defined by this vector. Mean-square fluctuations, 
J- 2 (0) ~ q 2 P 2 /( 2 n), is then calculated, where 9 is 2iv/q con¬ 
verted into arcsec. Mean-square fluctuations are divided by 
the fraction of unsubtracted pi xels in the ima g e to compen¬ 
sate for area lost by masking dArendt et al.ll201Cl ). Unless 
otherwise noted, error bars shown in figures represent one 
standard deviation of the Poisson errors, P 2 /\/N q , where 
N q is the number of Fourier elements in the associated bin. 
Because fluctuations may be reported as the mean-square 


2 http://archive.stsci.edu/pub/hlsp/udf/nicmos- 
treasury / table/szalayf 160wnfdet .fit 

3 http://archive.stsci.edu/pub/hlsp/udf/nicmos- 
treasury / table / newudfcommatable.txt 
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Original Mask 
Expansion Factor of 1.5 
Expansion Factor of 2.0 
■ Expansion Factor of 3.0 
Attenuated Border 


Figure 2. Example of mask expansion of the source outlined in Fig[l] (a): Original source subtraction showing significant residual wings, 
(b): The mask area is enlarged by a factor of three, (c): Source masks for expansions of 1.5, 2.0 and 3.0. Masked pixels are assigned 
values of zero in Panels (a) and (b). 



Figure 3. F160W images used to estimate effects of expand¬ 
ing individual source mask areas, (a): Original, source-subtracted 
image before mask expansion, (b): Mask-expanded image with 
source mask areas expanded by a factor of three. When compared 
to Panel (a), residual circumferential flux is markedly reduced, 
(c): Control Image. Random source masks are placed in random 
locations to exactly match subtracted pixel count of the mask- 
expanded image shown in Panel (b). (d): Noise image. Unsub¬ 
tracted (non-zero) pixels in the control image are replaced with 
Gaussian noise. Circles outline three of many randomly placed 
source masks generated for control images (lower panels) that are 
not present in the original or mask-expanded images (upper pan¬ 
els). The upper circle shows that wings of identified sources are 
not masked in the control images. 

(jKashlinskv et aD 120121 ; I Fernandez et al.l l2Q12h or root- 
mean-square (rm s), \/g 2 P 2 /(27t) (fKashlinsky et alj [2005i; 
Thompson et al.l l2007al : iHeleason. Ricotti. fc Kashlinskvl 
2012), both of these are calculated when appropriate. 


3 FLUX FROM EXTENDED WINGS 

The contribution of unsubtracted flux in galaxy wings to 
total fluctuations is estimated by incrementally expanding 
subtracted regions (masks) for individual sources. While this 
expansion effectively eliminates residual flux in wings, the 
masking process itself can distort the Fourier power spec¬ 
tra. To account for this distortion, comparisons are made to 
control images with the same number of unsubtracted pixels 
(see Section E2J). 


3.1 Expansion of individual source masks 

To quantify contributions from galaxy wings over a range 
of brightness, all source masks are uniformly expanded in a 
predefined manner. This is possible because the extraction 
process identifies specific objects, thereby allowing mask ex¬ 
pansion tailored to individual source size. Source mask ex¬ 
pansions are evaluated over a range of 1.1 to 6.0 times the 
original source-mask pixel count. Fig. 0 demonstrates this 
process applied to the single galaxy outlined in the bot¬ 
tom panels of Fig. |T] In the left panel, residual flux is seen 
extending beyond the original subtracted region. The mid¬ 
dle panel shows that residual flux from the same object is 
markedly reduced after the mask area is tripled. The right 
panel demonstrates the expansion process in more detail for 
expansion factors of 1.5, 2.0 and 3.0. Individual source masks 
are circumferentially extended in one pixel increments until 
the number of masked pixels is the desired multiple of its 
original value. This process will significantly over-subtract 
small objects where one pixel represents more than 10 per 
cent of the desired mask area. This problem is minimised by 
expanding the outermost portion of the mask on an image 
resampled with pixels 1/10 their original edge length using a 
cubic spline interpolation. These outermost pixels are then 
rebinned back to their original size by averaging over the 
100 resampled pixels. As shown in the right panel of Fig. [2] 
this results in an attenuation of these outer pixels rather 
than them being set to zero. The effect of this attenuation 
is included when calculating the total number of masked pix- 
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2n/q (arcsec) 2n/q (arcsec) 2n/q (arcsec) 

Figure 4. Mean-square fluctuations for unexpanded, control, and mask-expanded F160W images after expanding source-mask areas by 
a factor of four. Gaussian noise fluctuations are estimated separately for each image. Error bars represent lcr of Poisson errors. Shaded 
areas correspond to fluctuations above noise. 


Table 1. Percentage of masked image pixels for selected expan¬ 
sions. Values include effects of attenuated pixels (see Section l3.ll . 


Expansion factor 

Per cent of masked pixels (%) 

1.0 (baseline) 

13 

2.0 

22 

3.0 

30 

4.0 

38 

5.0 

44 

6.0 

50 


els on an image. The equivalent percentage of image pixels 
masked for a range of expansion factors is shown in Table [T] 
An example of an F160W image where source mask sizes 
are expanded by a factor of three is shown in Fig. (3}z>. When 
compared to the baseline source-subtracted image shown in 
Fig- Gib it is apparent that much of the residual flux in wings 
is eliminated. 

3.2 Generation of control images 

The masking process itself can cause a distortion of fluctua¬ 
tion s pectra that must be considered when making compar¬ 
isons (lArendt et all l20ld : iKashlinskv et all 120121 . Images 
are masked by multiplying the image field by a separate 
mask with pixel values of one everywhere except for source- 
subtracted regions that have pixel values of zero or, as shown 
in Fig. 2, a transmission factor. The Fourier transform of 
this multiplication is equivalent to the convolution of the 
Fourier transform of the mask with that of the unmasked 
image. Spectral properties of the mask are influenced by 
both the locations and sizes of individual source masks. The 
larger, brighter sources are not evenly distributed through¬ 
out the NUDF field and tend to be grouped towards the left 
side of the image (Fig. 1). Masking this uneven population 
of bright sources would be expected to generate some re¬ 
distribution of fluctuation power that is not representative 
of the true background spectrum. Because of these effects, 
mask-expanded spectra are compared to those of control im¬ 
ages generated by adding randomly placed source masks 
to the original source-subtracted NUDF image. This ad¬ 


ditional masking is done by randomly selecting individual 
source identification numbers from among the 4276 objects 
identified in the source extraction process. The mask from 
the selected object is then placed in a random location on 
the image. This process continues until the total number of 
masked pixels in the control image exactly matches that of 
its associated mask-expanded image. During this process it 
is important that the randomly placed source-masks remain 
on background pixels and do not subtract galaxy wings. This 
is accomplished by not allowing a randomly placed mask to 
encroach on a region encompassing four times the original 
mask size of an identified source. A region this size encloses 
essentially all of the residual wing flux associated with an 
object (see Section G-flll . An image created in this manner is 
shown in Fig. [3}i. Because locations of the randomly placed 
source-masks will affect fluctuation calculations, 1000 real¬ 
isations of the control images are generated. Contributions 
of noise to fluctuation power are estimated by replacing all 
nonzero pixels of an image with Gaussian noise of the same 
mean and standard deviation as background sky (Fig. (3jl). 
This is done for all images. Averaged fluctuation spectra of 
the 1000 control images and respective noise images are used 
for further calculations. 

3.3 Effects of mask expansion on fluctuations 

Figs [4] and [5] show results from 1.6 //m images after ex¬ 
panding source-mask areas by a factor of four. Mean-square 
fluctuations and associated Gaussian noise in unexpanded, 
control, and mask-expanded images are plotted in Fig. [4] 
Net fluctuations, calculated as 

J r L(e) = J r2 (e)-J r Lue) (i) 

are also shown for each of these images, where symbols 
are defined in Table [2] When compared to spectra from 
the unexpanded and control images, the image with source- 
expanded masks (right panel) has decreased fluctuations, 
which are most apparent above ~3 arcsec. Direct compar¬ 
isons of total mean-square fluctuations for these images are 
shown in Fig. 0 As shown in the upper panel of Fig. [|J] pro¬ 
gressive mask expansion causes mean-square fluctuations to 
noticeably decrease from about 1-50 arcsec relative to the 
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2n/q (arcsec) 

Figure 5. Mean-square fluctuations for unexpanded, control, and 
mask-expanded F160W images after expanding source-mask ar¬ 
eas by a factor of four. Error bars represent 1 rj of Poisson errors 
in the control image. 




Figure 6. Mean-square fluctuation changes at 1.6 fJ .m relative to 
unexpanded, source-subtracted image for four different expansion 
factors. Changes for source-expanded images are shown in the 
upper panel and those for control images in the lower panel. 


unexpanded image with a maximum reduction of almost 60 
per cent in the 10 to 40 arcsec range at an expansion fac¬ 
tor of six. While source-mask expansion generally results in 
a decrease in fluctuations, Figs [5] and [6] show an increase 
at ~80 arcsec where fluctuations in the mask-expanded im¬ 
ages exceed those of the unexpanded image. Although the 
actual difference is small, and on the order of Poisson noise, 
its effect on a percentage basis is exaggerated by the low 
fluctuation power at large scales. 

In theory, control images generated by randomly placing 
individual source-masks on the NUDF background should 
not alter the underlying power spectrum. However, the bot¬ 
tom panel of Fig. [5] shows that even after correcting for the 
fraction of masked pixels, the randomly placed source-masks 
on control images cause distortions in the spectral power 



1 10 100 

2-7r/q (arcsec) 


Figure 7. Mean-square and rms fluctuation changes at 1.6 (im 
relative to control images for four different source mask expansion 
factors. Inner error bars represent one standard deviation for 1000 
calculations for an expansion factor of four. Outer error bars also 
include Poisson errors. 

which is primarily due to the masking process itself (see 
Section [3721) . Per cent changes become progressively more 
significant at larger angular scales as the expansion factor is 
increased and can be as high as almost 40 per cent for an 
expansion factor of six. 

3.4 Contributions of unsubtracted galaxy wings 
to fluctuation power 

Source-mask expansion multiples of 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 
1.7, 1.8, 1.9, 2.0, 2.5, 3.0, 4.0, 5.0 and 6.0 in both filters are 
tested. At the maximum mask expansion of 6.0, 50 per cent 
of image pixels are set to zero. To compensate for spectral 
distortions caused by masking, particularly at high expan¬ 
sion factors, fluctuations in mask-expanded images are com¬ 
pared to the mean of their 1000 associated control images. 
Fig. [7] shows relative changes after expanding source-mask 
areas in the 1.6 /im image. The per cent change in mean- 
square fluctuation power at angle 9, A p(9), is calculated 
as 

A P(0) = 100 X ~ -^control (g) , (2) 

''control v / 

and for rms fluctuations as 

A p(6) = 100 x • F " xp( ^ ~ . (3) 

Relative changes plateau with increasing source mask ex¬ 
pansion and essentially no benefit is seen after an expansion 
of 4.0. The peak effect at 1.6 /im is at ~10 arcsec with de¬ 
creases of 52 ± 6 per cent for mean-square and 31 ± 3 per 
cent for rms fluctuations. The maximum changes at 1.1 /im 
are at 40 arcsec with decreases of 50 ± 25 per cent for mean- 
square fluctuations and 30 ± 11 per cent for rms values. 
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Table 2. Parameter Definitions 


Abbreviation Definition 


F sx P (0) 

F control (0) 

HO) 

F noise (,0) 

Fneb(0) 


Fluctuations in the image with expanded source masks 
Fluctuations in the control image. See Section 13.21 
Total fluctuations in an image, including noise 

Fluctuations in a ’’noise image” with unmasked pixels in original image replaced by Gaussian noise 
Difference between fluctuations in the original image and its associated noise image 


All fluctuations are a function of the angular scale, 0 



Figure 8. Integrated fluctuation changes as a function of source 
mask expansion for F160W and F110W images. Error bars in¬ 
clude both one standard deviation for the 1000 calculations and 
Poisson errors. 


Fig. [8] shows results for all expansions after summing 
fluctuations above 0.5 arcsec. The per cent decrease in mean- 
square fluctuations, A pi, at expansion factor, z, is 


Api = 100 x 


E^ 2 con tr ol(0j)-£^ 2 exp(6>j) 

j_j_ 

EE 2 control ( 0 j) 

j 


(4) 


and for rms fluctuations, 


Api = 100 X 


E K control W) - JEE 2 exp (*j) 


EE 2 control (*j) 


(5) 


As suggested in Fig. [TJi, subtracting the F110W wings has 
less effect than subtracting those of the F160W images. No 
significant benefit is noted in F160W images after individual 
mask areas are expanded by a factor of ~4.0 with only minor 
changes above 2.5. Eliminating residual flux in the wings of 
identified objects reduces total summed mean-square fluc¬ 
tuation power in 1.6 pm images by only 6.5 ± 0.3 per cent 
and rms fluctuations by 3.6 ± 0.2 per cent. Changes are sig¬ 


nificantly less at 1.1 pm with reductions of 1.4 ± 0.3 and 
0.7 ±0.1 per cent, respectively. 

As noted above, subtracting residual flux in the 
1.6 pm images reduces mean-square and rms fluctua¬ 
tions by as much as 52 and 31 per cent (Fig. 0. 
The significance of these decreases is estimated for pub¬ 
lications of four theoretical models of NIRB fluctua¬ 
tions in which the authors have cited NUDF results to 


support their conclusions (Heleason, Ricotti. & Kashlinskv 

2012); 

Fernandez et al. 20121; lYue et al. 2013al; Coorav et al. 

2012a 

). Fig. [9] compares fluctuations at 1.6 fin 1 after a 


source-mask expansion of 4.0 to the original NUDF data 
used in these publications. Data points from cited plots 
were extracted using online digitising software 0 While an 
expansion of 4.0 effectively eliminates unsubtracted flux in 
galaxy wings (Fig. 0, 38 per cent of total pixels are set 
to zero. As discussed in Section m and shown in Figs [6] 
and 0 if spectral distortions caused by the masking pro¬ 
cess itself are ignored, changes will be overestimated at 
angular scales larger than ~10 arcsec. Comparing fluctu¬ 
ation spectra of source-mask expanded images to control 
images with the same number of subtracted pixels provides 
a better compensation for masking effects than only cor¬ 
recting for source-subtracted areas. Therefore, mean-square 
and rms fluctuations for the mask-expanded images are cal¬ 
culated by multiplying their original, unexpanded values by 
F e 2 xp (6»)/F c 2 ontrol (0) and F e xp(0)/.^control(#), respectively. As 
shown in Fig. 0 subtracting residual flux from source wings 
would not affect any published conclusions drawn from the 
previously reported NUDF fluctuations. 


4 DISCUSSION AND CONCLUSIONS 

This study has shown that faint wings of source-subtracted 
galaxies account for only a small portion of total NUDF 
background fluctuation power above 0.5 arcsec. These re¬ 
sults are generally consistent with estimates of the con¬ 
tributions of residual galaxy wi ngs to spectral power 
in Spitzer /IRAC deep images IIKashlinskv et ahl 120051 : 
lArendt et aliFioTol f . These prior studies did not identify in¬ 
dividual sources and galaxy wings were removed by eq uall y 
increasing the size of all masked regions (lArendt et all 2010l. 
Radial mask expansion was limited to two pixels in order to 
retain enough unmasked pixels for reliable Fourier analy¬ 
sis. The higher resolution of the NUDF observations used in 
the present study allows residual flux to be eliminated using 


4 http://arohatgi.info/WebPlotDigitizer 
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27 i/q (arcsec) 




Figure 9. Fluctuations after subtracting source wings compared to the original (unexpanded) NUDF data for four published theoretical 
mo dels of the 1.6 Atm NIRB. (a): The s haded region corresponds to t he bracketed high (HFE ) and low (LFE) faint end areas shown in fig. 
9 of lHelgason, Ricotti, fc~K ashlinskv (2012). (b): Simulations from Fernandez et al.1 (120121 1 adapted from their fig. 4, include minimum 
masses of 10 s Mg (solid line) and 10® Mg (dotted line) without suppression as well as models with sm all haloes suppressed at high 
(dashed lines) and low (dash-dotted lines) efficie ncies, (c): Models shown from fig. 5 oflYue et all ll2013al l include contributions to the 
NIRB of their estimates for z > 5 (dashed lines!. iHel gason. Ri cotti, &: Kashlinskvl ll2012ll est imates for z < 5 (das h-dotted line), and the 
calculated total (solid line), (d): The models of halo occupation distributions from fig. 12 of lCoorav et al.l ll2012al l represent a threshold 
mass of 10 6 Mg with efficiency fractions of 0.03, 0.02, and 0.04 and associated ionisation photon escape fractions of 0.5, 0.9, and 0.1. 
These are shown as the solid, dashed, and dotted lines, respectively. Plot dimension ratios, axis labels, and plot symbols are similar to 
those in the original publications. Error bars in the NUDF data represent Poisson noise in the original (unexpanded) images and in many 
cases are obscured within the plot symbols. 


an algorithm that expands individual source masks propor¬ 
tional to the size of the underlying galaxy. This reduces the 
number of expanded pixels and allows a more robust analy- 
sis of the effects o f gala xy wings on measured fluctuations. 
iKashlinskv et ah (2005 1 found a decrease of less than a few 
per cent in Pilq ) after r adially expanding masks by two 
pixels while lArendt et al.l (|2010| ) noted that in most cases 
the largest changes in rms fluctuations were less than the 
lcr uncertainties. Both of these findings are significantly less 
than the decreases observed in this study (Figs [7] and [9} and 


probably result from a more complete source-subtraction in 
their original images. 

A limitation to the methodology used in this study is 
the degree to which control images can compensate for spec¬ 
tral distortions caused by source-mask expansion. Individ¬ 
ual source-masks in control images have a smaller average 
size and a more even distribution than those in the source- 
mask expanded images and, therefore, will affect underly¬ 
ing fluctuations differently. Although using control masks 
as generated in this study will not fully correct for spec¬ 
tral distortions caused by mask expansion, they do provide 
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a better compensation than only dividing by the fraction 
of unmasked pixels. This is important because, as shown in 
Figs [h] (upper panel) and[7| both the magnitude of changes 
and the expansion factor required to eliminate wings would 
be overestimated if control images were not used. 

The IHL model for fluctuations applied to observations 
from both the Spitzer Deep, Wide-Field Survey (Cooray et 
al. 2012b) and CIBER (Zemcov et al. 2014) predicts effects 
becoming significant at angular scales larger than those en¬ 
compassing the NUDF field. Yue et al. (2013b) fit their 
DCBH model for fluctuations at 3.6 and 4.5 /rm to scales 
> 100 arcsec and, as suggested in fig. 4 of their paper, 
DCBHs would have negligible contributions at 1.6 and 1.1 
/im. Therefore, studies postulating IHL or DCBHs to ac¬ 
count for the NIRB did not use NUDF results to support 
their hypotheses. 

As shown in Fig. [9] even the maximum decreases 
for mean-square and rms fluctuations are barely notice¬ 
able on a log scale. Results after eliminating source wings 
from NUDF images represent the upper limits of fluctu¬ 
ations at 1.1 and 1.6 fim and continue to remain well 
above th ose predicted by models for known galaxy pop¬ 
ulations jHelgaso mjtkxrtti 1 _fcj£ashlmsk'\j_ 2012|) o r high- 
redshift galaxies ([Fernandez et al.l 120121 ; lYue et al.1 l2013al : 
ICoorav et all I2012al l. The revised fluctuation spectra pre¬ 
sented in this paper would not affect any of these conclu¬ 
sions. These results exclude previous suggestions that ex¬ 
tended wings of detected galaxies may be a major contrib¬ 
utor to fluctuations in source-subtracted NUDF images and 
confirm that other explanations are required to account for 
most of the NIRB fluctuations at the angular scales of NUDF 
observations. 
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